Field Induced Long-lived Super-Molecules 
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We demonstrate that the long-lived bound states (super-molecules) can exist in the dilute limit 
when we tune the shape of effective potential between polar molecules by an external microwave 
field. Binding energies, average sizes, and phase diagrams for both s-orbital (bosons) and p-orbital 
(fermions) dimers are studied, together with bosonic trimer states. We explicitly show that the non- 
adiabatic transition rate can be easily tuned small for such ground state super-molecules, so that 
the system can be stable from collapse even near the associated potential resonance. Our results, 
therefore, suggest a feasible cold molecule system to investigate both novel few-body and many-body 
physics (for example, the p-wave BCS-BEC crossover for fermions and the paired condensate for 
bosons) that can not be easily accessed in single species atomic gases. 



Since the successful realization of a high-phase-space 
density of molecules at JILA [lj, systems of ultracold 
polar molecules become one of the most exciting devel- 
opments in strongly correlated physics. However, the 
strong dipolar interaction between polar molecules may 
render the system unstable to collapse in the high density 
regime (also observed in atomic gases with negative scat- 
tering length). Recently there have been several methods 
proposed to overcome such an unwanted situation by re- 
ducing the system dimensions in an optical lattice 043- 
On the other hand, an effective repulsive core can also 
be generated to stabilize the system if a proper detuned 
AC field is applied [||. These mechanism open new pos- 
sibilities to investigate interesting many-body physics of 
polar molecules, including the p-wave superfluid between 
fermionic molecules the itinerant ferro-electricity 
and dimerized superfluidity [lOJ etc. 

In addition to many-body physics, it is also interest- 
ing to study the few-body physics of polar molecules in 
the dilute limit. In Ref. [11(, the authors proposed a 
DC-field linked bound state by pumping the electrons 
to a meta-stable state, which, however, has a very short 
life time and is therefore not easily accessible. DC field 
induced inter-layer (-tube) bound states in a multi- layer (- 
tube) system were also proposed Q, but the mutual 
long-ranged repulsion makes the critical density for such 
bound state formation much lower than the density avail- 
able in current experiments jl2j. 

In this paper, we investigate a different few-body 
bound state (super- molecule) in the presence of an exter- 
nal AC field [13]. The effective potential between polar 
molecules can be tuned to have an attractive and shallow 
well (see Fig. QJa) an d (b)), which supports only a few 
or even no bound state (dimers and trimers) between 
molecules. Their average sizes are of sub-micrometer 
scale, strongly suppressing the few-body inelastic scat- 
tering rate and stabilizing the whole system from col- 
lapse. As a result, our work suggests a new area for novel 
few-body physics with a long-lived ground state. Some 




FIG. 1: (Color online) (a) Rotational eigenstates of a single 
molecule. The microwave field couples states |0,0) and 1, 1) 
with the Rabi frequency SI and the detuning A (circularly 
polarized in the x — y plane as shown in the inset of (b)). 
(b) The effective interaction (14ff(r)) for different inclination 
angles 9. 9 = 0, tt/6, it/3, and 7r/2 respectively from top to 
bottom, (c) and (d) show the 3D contour plots of the ground 
state wavefunction of a dimer in the s- and the p^-state re- 
spectively. Red(blue) color stands for positive(negative) value 
of the wavefunction. Here we use Q,/ '£?=(]. 01 and A/B=0.01 
for SrO molecules. 



important many-body physics not accessible in atomic 
gases, for example the p-wave BCS-BEC crossover be- 
tween single species fermions [1J] and the Z2 transition 
[HI between the condensate and the pair condensate of 
single species bosons near the Feshbach resonance, may 
become feasible without particle loss. 

In the low temperature and dilute limit, the only 
important internal degrees of freedom is the rotational 
states, if considering close-shell molecules (SrO as an 
example) without hyperfine structure. As a result, 
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the general two-body Hamiltonian with dipolar interac- 
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tion becomes [8j|: if = Yli=i 2 

di-da-aCdx-eiaK^-e^) wh d £ are ^ molecu . 

lar mass, angular momentum operator and the rotational 
constant respectively. E(i) is the external field, and d; 
is the electric dipole moment operator, and Pi are the 
position and momentum operators. e±2 = ri2/|ri2| and 
ri2 = ri — Y2- Here we will consider only the case when a 
circularly-polarized microwave field is tuned to be close 
to the transition energy between |0, 0) and |1, 1) (see Fig. 
HJa)). Here |J, M) is the rotational state labeled by the 
angular momentum quantum number. 

Using the adiabatic and rotating wave approximations 
Q, we obtain the full effective interaction Q, which has 
the leading order terms in the long distance: 
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with d being the bare electric 
3 cos 2 9) 2 
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and Cf\e) = 
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with £ = ^/l + 4(f2/A) 2 . 9 is the angle between the 
microwave propagation and the relative position of the 
two molecules (see the inset of Fig. [ljb)) . We let h = 1. 
We can see that V e g(r) is dominated by the dipolar-like 
term with the effective dipole moment, d c e, as |r| — > 00, 
but its next order correction is like a hard-core potential 
to prevent possible chemical reactions (HQ. Therefore, 
one can always expect a potential minimum at |r| = r m 
in the x — y plane. 

In Fig. [lib), we show the full numerical V e g(r) for 
different 9 with the detuning A/B = 10~ 4 and the Rabi 
frequency VL/B = 1CT 3 . Note that, by using r B = (§-) 1 ^ 3 
and B as the unit of length and energy, Kff(r) is universal 
for all molecules, once the external field parameters (£l/B 
and A/B) are specified. Using SrO as an example (d = 
8.9 Debye B ~ 10 10 Hz), we find tb ~ Ifnm and hence 
the super-molecules can be as large as 100 nm as shown 
in Fig. [ljc) and (d). The dimer super- molecule has a 
donut-type wave function for bosons (s-orbital state) and 
has two knots for fcrmions in the p^-orbital state [13]. 

Since all the important and non-trivial properties of 
the effective potential lies on the x — y plane, from now 
on, we will consider 2D molecular motion for simplic- 
ity and use = U c ff(|r|,# = 7r/2) as the effective 
interaction. Such system can be physically realized by 
adding a strong confinement in the z direction. The fi- 
nite width correction can be neglected because the donut 
type super-molecule has a very small thickness (less than 
8vb ~ 88 nm even for SrO). 

In Fig. [UJa) and (b), we show the phase diagrams for 
the number of 2D dimer states in terms of the detuning 
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FIG. 2: (Color online) (a) and (b) are phase diagrams for the 
number of s- and p-orbital states of a dimer super-molecule. 

(c) is the resonance of the s-wave scattering length (0,3, see 
text) as a function of Q (i.e. the yellow vertical line in (a)). 

(d) is the corresponding binding energies for the s-orbital 
ground (blue solid line) and the first excited (orange solid 
line) states. Results for the p-orbital ground state, the bosonic 
trimer ground, the first excited and the second excited state 
are shown in green dashed line, blue circles, orange squares 
and purple triangles respectively. 



and Rabi frequencies, including both s- and p-orbitals. 
In the limit of small detuning (say A/B < I0~ 4 ), the 
bound state is solely determined by f2, while in the limit 
of large detuning, it depends on A only. In the inter- 
mediate regime, there is a triangular window (0.001 < 
A/B < 0.01), where no dimer bound state exists in both 
orbital. 

The appearance of a two-body bound state implies a 
potential resonance in the low energy scattering ampli- 
tude. For our 2D effective potential with a r~ 3 tail in 
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where ip(r) is the relative wavefunction, /3q = e 7 
and 7 ~ 0.57722 is Euler's constant. Here 03 = 
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is the effective scattering length [if 



with Pq being the only parameter determined from the 
short-range effective interaction (see Refs. [16| for de- 
tails) . The incident momentum k is related to the system 
density. The 2D pseudo-potential is valid for ka 3 ^ 1. 

In Fig. djc), we show the calculated s-wave scatter- 
ing length, 03, as a function of the Rabi frequency at 
A/B = 0.003. The corresponding binding energies are 
also shown in Fig. [UJd). As expected, the divergence of 
the scattering length (similar to Feshbach resonance) in- 
dicates the appearance of an s-orbital bound state with 
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a binding energy typically smaller than one [iK (taking 
B ~ 10 GHz). The tunability of the scattering length 
will certainly affect the many-body properties, and this 
result remains true for the 3D system. 

In addition to dimer states, we also investigate three- 
body bound states (trimers). We note that the Efimov 
related physics in 3D has recently been investigated for 
dipolar atoms [I?} or 2D systems [l8| with short-ranged 
interaction. However, the effective interactions of our 
system are not simple dipolar-like, but with several in- 
trinsic length scales. Therefore one should not expect 
universal results in our current system. Here we apply 
the stochastic variational method (SVM) with Gaussian 
function basis to study trimer states of bosonic molecules 
pj| . The uncertainty of our calculation is about three or- 
ders of magnitude smaller than the binding energies. In 
Fig. [3Ja) and (b), we show the typical trimer wavefunc- 
tion and its size as a function of the Rabi frequency. 

In Fig. [2jd), we also show the binding energy of a 
typical trimer super-molecule as a function of Rabi fre- 
quency. If we stay close to the threshold of bound state 
formation, there is a finite parameter range (2 x 10~ 4 < 
Q/B < 4 x 10~ 4 ) in which the energy differences from 
dimer's are less than 10 _6 i? ~ lOfc Hz for SrO, smaller 
than the typical trapping frequency. In other words, 
even though two super-molecules (four particles) may 
have inelastic scattering, leading to the formation of one 
trimer and one free particle, the energy gain of the later 
should not be high enough to escape the trap. As a re- 
sult, we do not expect a significant loss rate from such 
four-body (two super-molecule) scattering, and the re- 
sulting ground state should be a mixture of dimer and 
trimer super-molecules. This is a novel few-body system 
composed of field-dressed polar molecules. For molecules 
with smaller dipole moment, tb is smaller. The stable 
regime of mixture would be enlarged. When we move 
away from the threshold, the excited states start to show 
up (at fl/B - 5 x 1CT 4 with A/B = 0.003) with energy 
lower than a dimer. 

Now we consider the loss rate of a super-molecule due 
to the breakdown of the adiabatic approximation. The ef- 
fective potential is made of the avoided crossing between 
two eigenstates in the weak field limit (see Fig. UK a)) 
[1GJ]. They can be denoted to be |1) (effective potential) 
and |2), merge to \g, g, N) and \g, e+, N — 1) respectively 
as r — > oo. Here \g,g,N) denotes the two molecules in 
the rotational ground state with the photon number N, 
while \g, e+; N — 1) denotes the symmetric superposition 
of one molecule in the ground state and the other one in 
the excited state with one photon absorption. The non- 
adiabatic loss thus mainly from the inter-state transition 
from |1) to 1 2) due to the relative kinetic energy between 
two particles neglected before (see Fig. Ufa)). The transi- 
tion to other states is certainly much smaller and there- 
fore neglected here. Using Landau- Zener formula [20| . 
we can estimate the transition rate to be 72 = /aveC -2 
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FIG. 3: (Color online) (a) The wavefunction of a trimer state 
with A/B = 0.003 and Q/B = 7 x 10~ 4 with the coordinates 
defined in the inset: ri2 is the distance between two particles, 
and s is the distance between the third and the center of mass 
of the first two. Their relative angle has been integrated out. 
(b) The sizes of bosonic trimer (red circles), s-orbital and p- 
orbital dimers (blue solid and green dashed lines respectively) 
as a function of fl for A / B = 0.003. The vertical line indicates 
the appearance of the p-orbital dimer. 



with the transition parameter, 

r = 



\dE 12 (r)/dt\ 7 
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Here / avo is the frequency for two molecules to collide 
with each other, E 12 = E\\\ — E\%\ is the energy dif- 
ference between the two unperturbed energy eigenstates, 
and V12 — £1 is the coupling between them via microwave 
absorption. Within the semi-classical approximation, we 
have dEi2(r)/dt = E' 12 (r)v(r), where v(r) = dr/dt = 
y / 2(Eb — V2D (?■))//•* is the local velocity. Here Eb is the 
bound state energy and /1 = m/2 is the reduced mass. By 
estimating f~^ e = Jjj 2 for the two turning points at 

T\ and r 2 , we can calculate the lifetime t 2 = j 2 1 directly. 
In Fig. Sfb) we show the full calculation of the transition 
parameter, T, as a function of fl for SrO (B ~ 10 10 Hz 
and tb ~ llnm). The typical value of T is always larger 
than 5 for SrO, and therefore the obtained life time (t 2 ) 
is very long. We also get t 2 ~ 0.1 sec for KRb molecule 
(B ~ 10 9 Hz and re ~ 3.5nm, if omitting its fine struc- 
ture). 

In order to understand the life time in more details, we 
can qualitatively estimate its order of magnitude when 
the bound state is just formed (Eb — > 0). First we have 
/avo ~ Sv/2Sr ~ m(6v) 2 /2 ~ Ek, where Sr(Sv ~ 1/mSr) 
is the uncertainty of radial position(velocity) and Ek ~ 
V m is the kinetic energy. Since £12 (r) — > A as r — > 00 
and changes within the length scale of r rn ~ 5r/C with a 
field-dependent constant C, we then have E' 12 (r) ~ A/r m 
and therefore T ~ C |fi| 2 /AEk ■ For example, for the 
parameter of SrO with fl/B ~ 10~ 4 and A/B ~ 10~ 3 , 
we have C ~ 1 (see Fig. Ufa)) and hence a very long 
life time as calculated above (Ek/B ~ 4 x 10~ 7 and 
r ~ 23). We note that, different from the situation in 
the continuous state, the life time of a super-molecule 
is longer for stronger dipole moment because the length 
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FIG. 4: (Color online) (a) shows the non-adiabatic transition 
of a dimer state. The thick blue(thin red) lines are eigenstates 
|1)(|2)). Dashed lines are states as — > 0. (b) The calculated 
transition parameter, F, for the s-orbital (solid line) and p- 
orbital (dashed line) dimer ground state. A/B = 0.001, 0.002, 
and 0.003 respectively from top to bottom. The vertical line 
indicates the appearance of the p-orbital dimer for A/B — 
0.003. 

scale, tb = {d 2 /B) 1 ^, as well as the radial size of the 
super-molecules, also becomes larger, reducing the rel- 
ative collision velocity between two molecules. This is 
consistent with the fact that the adiabatic approximation 
is more accurate for particles with larger mass. Similar 
analysis can be also applied to the trimer case. 

The results presented here open a new area for study- 
ing few-body physics as well as many-body physics with 
ultracold molecules. In usual atomic gases, Feshbach 
molecule is the highest energy bound state with many 
lower energy states far below. As a result, the whole sys- 



tem of the same species is usually unstable to collapse 
due to the three-body inelastic collision in the s-wave 
(p-wave) channel for identical bosons (fermions). This 
fact makes the experimental realization of a Zi quantum 
phase transition [1 51 ] between a condensate and a pair 
condensate of bosons near a Feshbach resonance almost 
impossible. Similar conclusion also applies to the p-wave 
crossover, or the p z -to-p x + ip y transition, predicted in 
polarized fermion systems [14j . In the polar molecule 
system we discuss here, however, the ground state super- 
molecules is long-lived almost without two- and three- 
body loss near the threshold of dimer formation, sug- 
gesting a much more feasible scenario to study important 
many-body physics such as the phases mentioned above. 
Although larger clusters may exist when the density is 
higher, above statement and prediction should still ap- 
ply near the threshold of dimer-trimer formation. Much 
more complicated analysis will be needed to fully address 
this issue. 

In summary, we show that the shapes of the 2D effec- 
tive potentials between polar molecules can be designed 
to make long-lived super-molecules. This creates a fea- 
sible system for studying novel few-body physics as well 
as many-body problems. 
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